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We use three-dimensional molecular dynamics simulations of magnetically interacting pancake vortices to study 
the dynamic phases of vortex lattices in highly anisotropic materials such as BSCCO. Our model treats the 
magnetic interactions of the pancakes exactly, with long-range logarithmic interactions both within and between 
planes. The pancake vortices decouple at low drives and show two-dimensional plastic flow. The vortex lattice 
both recouples and reorders as the driving current is increased, eventually forming a recoupled crystalline-like 
state at high drives. We construct a phase diagram as a function of interlayer coupling and show the relationship 
between the recoupling transition and the single-layer reordering transitions. 



In highly anisotropic superconductors such as 
BSCCO, the vortex lattice is composed of individ- 
ual pancake vortices that may be either coupled 
or decoupled between layers depending on such 
factors as the material stoichiometry or the mag- 
nitude and angle of the applied magnetic field. 
Of particular interest is the possible relationship 
between a coupling/decoupling transition and the 
widely studied second peak or fishtail effect M . 

As a function of interlayer coupling strength s, 
there are two limits of vortex behavior in a system 
containing pointlike disorder. For zero interlayer 
coupling s = 0, each plane behaves as an indepen- 
dent two-dimensional (2D) system. For infinite 
interlayer coupling s = oo, the vortices form per- 
fectly straight three-dimensional (3D) lines, and 
all of the planes move in unison. At finite cou- 
pling strength ,s ^ 0, a transition between these 
types of behaviors should occur with coupling 
strength, but it is unclear whether this transition 
is sharp or if an intermediate state of the lattice 
exists. Furthermore, it is known that 2D systems 
with pointlikc pinning can exhibit dynamic re- 
ordering under the influence of an applied driving 
current, passing from a liquid- like state at zero 
drive to a recrystallized state at high current ||- 



|J. Thus, in a 3D system, a dynamically driven 
recoupling transition could be expected, but it is 
unclear where this transition falls in relation to 
the 2D reordering transitions already seen. 

To study the coupling transitions, we have de- 
veloped a simulation containing the correct mag- 
netic interactions between pancakes This 
interaction is long range both in and between 
planes, and is treated according to Ref. 

The overdamped equation of motion, T = 0, 
for vortex i is given by fj = Y2j=i VU(py, Zij) + 
f i p + fd — v^, where N v is the number of vor- 
tices, p and z are the distance between pancakes 
in cylindrical coordinates. The magnetic energy 
between pancakes is 
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where R = y/ z 2 + p 2 , Ei(x) — J*°° exp(~x/X)/p' 
and eo = $o/(47r£) 2 . The pointlike pins are ran- 
domly distributed in each layer and modeled by 
parabolic traps. We vary the relative strength of 
the interlayer coupling using the prefactor s. We 
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have simulated a 16A x 16A system containing 89 
vortices and 4 layers, with a total of 356 pancake 
vortices. Further work on systems containing up 
to 16 layers will be reported elsewhere 0. 

In Fig. 1(a) we present a phase diagram as 
a function of interlayer coupling strength s and 
driving force fd- At zero drive, we find a recou- 
pling transition for a coupling strength of s > 4.5. 
In samples with s > 5, the pancakes remain cou- 
pled into lines at all drives and show the same 
transitions seen in previous work |3), exhibiting 
plastic flow of stiff lines above depinning, and re- 
ordering into first a smectic state and then a re- 
crystallized state of stiff lines at higher drives. 

For samples with weaker interlayer coupling, 
s < 5, the vortex lattice is broken into decoupled 
planes at zero drive. Upon application of a driv- 
ing current, the samples exhibit 2D plastic flow in 
which each layer moves independently of the oth- 
ers. Once the individual layers reach the driving 
force at which a transition to a smectic state oc- 
curs, the vortices simultaneously form the smectic 
state and recouple, as can be seen from the mea- 
sure shown in Fig. 1(c) of the z-axis correlation 
C z = l-<(r i) £-r< )i+ i)e(ao/2-|(r i) £-r i ,i + i)|)) J 
where ao is the vortex lattice constant. The dy- 
namic recoupling transition line follows the smec- 
tic transition line down to s = 2 and is associated 
with a peak in the dV/dl curve seen in Fig. 1(b). 
Both the static and dynamic transition lines be- 
tween decoupled 2D and recoupled 3D behavior 
are sharp. 

As a function of the number of layers, we ob- 
serve the same behavior, but the depinning cur- 
rent in the 3D stiff state drops with the number of 
layers. The depinning current in the 2D decou- 
pled phase is not affected since in this case the 
individual planes behave as isolated entities. The 
recoupling transition sharpens as the number of 
layers is increased 0. 
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Figure 1. (a) Phase diagram for varying inter- 
layer coupling s and driving force fd- Circles: de- 
pinning line; diamonds: decoupling line; squares: 
smectic transition line; triangles: recrystalliza- 
tion line, (b) V x (circles) and dV/dl (+ signs) for 
s = 2.0. (c) C z for s = 2.0 (d) V x (circles) and 
dV/dl (+ signs) for s = 8.0. (e) C z for s = 8.0 
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